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ABSTRACT 

For quark matter studies in astrophysics the thermodynamic bag model (tdBAG) has been widely 
used. Despite its success it fails to account for various phenomena expected from Quantum-Chromo- 
Dynamics (QCD). We suggest a straightforward extension of tdBAG in order to take the dynamical 
breaking of chiral symmetry and the influence of vector interactions explicitly into account. As for 
tdBAG the model mimics confinement in a phenomenological approach. It is based on an analysis 
of the Nambu-Jona-Lasinio (NJL) model at finite density. Furthermore, we demonstrate how NJL 
and bag models in this regime follow from the more general and QCD based framework of Dyson- 
Schwinger (DS) equations in medium by assuming a simple gluon contact interaction. Based on 
our simple and novel model, we construct quark hadron hybrid equations of state (EoS) and study 
systematically chiral and deconfinement phase transitions, the appearance of s-quarks and the role of 
vector interaction. We further study these aspects for matter in /3-equilibrium at zero temperature, 
with particular focus on the current ~ 2 M 0 maximum mass constraint for neutron stars. Our 
approach indicates that the currently only theoretical evidence for the hypothesis of stable strange 
matter is an artifact of tdBAG and results from neglecting the dynamical breaking of chiral symmetry. 
Subject headings: dense matter — equation of state — elementary particles: quarks — stars: neutron 


1. INTRODUCTION 

QCD is believed to be the correct theory of strongly in¬ 
teracting matter. It inspires multifaceted research, both, 
theoretical and experimental (Brodsky et al. (2015) 
highlight this statement concerning the impact of QCD 
on hadron physics). Lattice QCD as the ab-initio ap¬ 
proach to solve QCD numerically is increasingly success¬ 
ful in vacuum, at finite temperatures and at small chem¬ 
ical potentials (cf. Fodor & Katz 2004; Aoki et al. 2006, 
and references therein). It is not suited to provide so¬ 
lutions for low temperatures at chemical potentials in 
vicinity of the predicted deconfinement phase transition. 
This is the domain of interest for astrophysical applica¬ 
tions, e.g., simulations of neutron stars (NS) and core¬ 
collapse supernovae. The influence of physically realized 
deconfined quark matter in NS/protoneutron stars and 
supernovae on potential observables is a long standing re¬ 
search topic (cf. Alford et al. 2007; Fischer et al. 2011; 
Klahn ct al. 2006, 2007, 2012; Nakazato et al. 2008; 
Pagliara et al. 2009; Pons ct al. 2001; Sagcrt ct al. 2009; 
Schertler et al. 2000). 

As an alternative to lattice calculations the DS ap¬ 
proach solves QCD’s gap equations within a given trun¬ 
cation scheme. It is applicable in the whole temperature- 
density domain of the QCD phase diagram. Al¬ 
though the DS formalism has been successfully applied 
to understand hadron physics on the quark level (cf. 
Bashir et al. 2012; Cloet & Roberts 2014; Chang et al. 
2011; Roberts 2012, and references therein) only little 
work has been done to exploit it at finite densities (cf. 
Roberts & Schmidt 2000) and to describe the EoS of de¬ 
confined quark matter (cf. Chen et al. 2008, 2011, 2015; 
Klahn et al. 2010). These explorative studies promise 
deeper insights into the EoS of cold, dense matter in a 
domain characterized by non-perturbative QCD. 


Perturbative QCD, valid near the limit of asymptotic 
freedom where quarks are no longer strongly coupled pro¬ 
vides a valuable benchmark for DS in-medium studies. 
Recent work by Kurkcla et al. (2014) illustrates how per- 
purbative QCD can be used in a reasonable way to put 
asymptotic high-density constraints on the EoS for com¬ 
pact stars. 

The majority of studies of dense quark matter in as¬ 
trophysical systems is based on effective models. These 
mimic certain but not necessarily all key features inher¬ 
ent to QCD, e.g., the breaking/restoration of chiral sym¬ 
metry and the effect of deconfinement on the EoS and 
related quantities. 

In this work we demonstrate that the two most com¬ 
monly used effective models in astrophysics, namely td¬ 
BAG as introduced in Farlii & Jaffe (1984), and mod¬ 
els of the NJL type (cf. Nambu & Jona-Lasinio 1961; 
Klevansky 1992; Buballa 2005) can both be understood 
as solutions of QCD’s in-medium DS gap equations 
within a particular set of approximations. From this per¬ 
spective we suggest modifications to tdBAG which make 
it consistent with the standard NJL approach: our vector 
interaction enhanced bag model (vBAG) accounts in a 
parameterized form for (a) the flavor dependent restora¬ 
tion of chiral symmetry, (b) modifications of the EoS 
due to vector interactions, and (c) - in the spirit of the 
tdBAG - a phenomenological correction to the EoS due 
to the deconfinement transition which one could other¬ 
wise not account for. Once parameterized, vBAG is as 
convenient to apply as the original tdBAG and therefore 
well suited for subsequent studies of dense and decon¬ 
fined quark matter at finite densities and temperatures 
in astrophysics. 

The manuscript is organized as follows: In sec. 2 we 
review briefly the in-medium DS gap equations which 
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we then reduce to the NJL model by assuming a simple 
contact interaction for the gluon propagator in sec. 3. In 
sec. 4 we begin to introduce our vBAG model and illus¬ 
trate the small impact of the residual dynamical chiral 
symmetry breaking (DySB) on the EoS in the chirally 
nearly restored phase in the NJL model. We also point 
to the importance of the flavor dependence of the chiral 
transition. In sec. 5 we add vector interactions to the 
model and obtain a set of equations which distinguish 
our vBAG from the standard tdBAG. We then discuss 
a setup for a phase transition construction which can 
ensure a simultaneous chiral and deconfinement phase 
transition in sec. 6, and apply the resulting quark-hadron 
hybrid EoS to describe high mass compact stars in sec. 7. 
In sec. 8 we review the long standing hypothesis of abso¬ 
lutely stable strange quark matter as the ground state of 
matter, and finally conclude with a summary in sec. 9. 

2. DYSON SCHWINGER FORMALISM WITH A CONTACT 
INTERACTION 


The key quantity of our analysis is the in-medium prop¬ 
agator of a single quark flavor (cf. Rusnak & Furnstahl 
1995; Roberts & Schmidt 2000), 

p) = tfpA(p; p) + i^ A p A C{p] p) + B(p- p) , (1) 


with P 4 = P 4 + ip, where p is the chemical potential. 
The gap solutions A, B , and C are obtained from the 
gap equation 


S 1 (p;p) = ijp + ij4P4 + m + '£(p-p), (2) 


s ^ = Ja y a2Dpcj ^ _ ^ ’ 

(3) 

where m is the current-quark mass, D pa (p) the gluon 
propagator and T“(p;g) the quark-gluon vertex. A rep¬ 
resents a regularisation mass scale which, in a realistic 
treatment, would be removed from the model by tak¬ 
ing the limit A —> oo. For the NJL model this pro¬ 
cedure fails and A usually represents a simple cut-off 
for momentum integrals. Following the approach from 
recent vacuum studies, cf. Gutierrez-Guerrero et al. 
(2010), we apply the rainbow truncation T“(p; q) = 4f7 CT 
and assume a contact interaction in momentum space, 
g 2 D p(J (p — q) = 5 pa ^rQ( A 2 — p 2 ). The Heaviside func¬ 
tion 0 provides a 3-momentum cutoff for all momenta 
p 2 > A 2 in order to regularize ultraviolet divergences 
inherent in Eq. (3). Different regularisation procedures 
are available 1 , we chose this hard cut-off scheme to match 
our model with tdBAG, i.e. describe quarks as a quasi 
ideal gas of Fermions, me is a gluon mass scale which in 
this model simply defines the coupling strength. These 
approximations are sufficient to solve Eq.(2). The A-gap 
has only a trivial, medium independent solution, A = 1. 
The remaining gaps take the following form, 


_ 16 N c f d\ B q 

P m + 9 ml A (2tt)4 cpA 2 + + B\ 


( 4 ) 


_~2 , 8N c [ P4q4C q 

P 4 P Pi+ 9m 2 G J A (2Try<pA 2 + qtC 2 q +Bi ' ( ° } 

1 In fact, the regularisation scheme does not have to affect ul¬ 
traviolet divergencies only. E.g., the very infrared, gluon mediated 
behavior is ill described in NJL type models; an IR cutoff scheme 
can remove unphysical implications (Ebert et al. 1996). 


The integrals do not explicitly depend on the external 
momentum p. Consequently, both gap solutions are con¬ 
stant and vary only with p. One easily verifies from 
Eq.(5) that the vector gap C induces a medium depen¬ 
dent shift u> of the chemical potential p, 

P 4 C = P 4 + iui = P 4 + ip* = P4 , (6) 


where the effective chemical potential p* is introduced 
in order to maintain a quasi particle description of the 
quarks. The vector condensate oj is given by 


_ * _ 8 N c r d 4 q _ iq 4 

9 m 2 G J A (2tt) 4 q 2 + qj + B 2 


( 7 ) 


Instead of solving this equation to determine p* for a 
given p it is more convenient to fix an arbitrary effective 
chemical potential p* and determine the actual chemical 
potential p post priori. Then only the remaining Eq. (4) , 
rewritten in terms of the effective chemical potential p *, 
needs to be solved self-consistently, 


„ :17 1 piew /• A 1 

9 m G J A (2tt) 4 q 2 + q 2 + B 2 


( 8 ) 


In analogy to Eq.(7) the scalar condensate is defined as 
cf> = B — m. As we aim to provide an EoS at finite densi¬ 
ties and temperatures we apply Matsubara-summations 
to obtain the following set of equations, 


B = m 


4W 

9 m G 


n s {p*,B) 


* 2iV c 

d ’=> 1 + 7ttM/ 4 i-Sb 
9 rrin 


( 9 ) 

( 10 ) 


with the single-flavor scalar and vector densities, n s and 
n v , of an ideal spin-degenerate Fermi gas, 


■=*£ 
± • 

■=*£ 


d 3 p B fl 1 > 

) ,(11) 

(2? r) 3 E \ 2 

1 + exp (E ± /T)) 

d 3 p 

Tl 

(12) 

(27r) 3 1 + exp (E ± /T) ’ 


with E 2 = p 2 + B 2 and E ± = E ± p*. 

It is common for DS calculations to determine the ther¬ 
modynamic pressure in steepest descent approximation. 
It consists of ideal Fermi-gas and interaction contribu¬ 
tions, 


Pfg =TrlnS' _1 

= 2A T c J ^ AL. \ n (f +p 2 + R2) ; ( 13 ) 

P 7 = ■-^TrES = ^m 2 G to 2 - |m. (14) 

The pressure of the free Fermi gas consists of vacuum 
and kinetic parts, Pfg = Pfg + Pfg 1 with 
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Fig. 1.— (color online) Single flavor dynamical masses (black) and corresponding pressure (red) computed within the NJL model. The 
latter is well fitted by the pressure of an ideal Fermi gas (with bare quark mass rrif ) shifted by a chiral bag constant B x (blue). Parameters: 
set IV, table 1. 


Divergences occur only in the vacuum contributions to Eqs. (20) exactly reproduce the gap Eqs. (10) derived 
the total pressure, given as within the DS framework, if one identifies 


Pvac = 2N c [ 

J a (2ti- r 

ml f r d 3 p B 0 \ 2 

9 m% \J A (2tt)3 E ) 


(17) 


where the 2 nd term on the right-hand side originates from 
the scalar condensate <j). Bq is the dressed quark mass at 
p = T = 0. The remaining medium dependent pressure 
is free of divergences and does, in a technical sense, not 
require regularisation. 


3. NJL MODEL COMPARISON 


We compare the previous results to a simple model of 
the NJL-type. The Lagrangian is separated into free and 
interaction parts, Cq +C\, where the latter accounts for 
scalar and vector interactions, 
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Ci = C s +Cy = Gs^^qfTaqf) 2 + G v (q f ij 0 qf) 2 , (18) 

a —0 



v — 



( 21 ) 


The relation of the NJL model coupling constants, G v = 
G s / 2, is consistent with the result obtained after Fierz 
transformation of the one-gluon exchange interaction (cf. 
Buballa 2005). We summarize the first part of this work 
by observing that the results we obtained within the 
DS approach under the assumption of a contact inter¬ 
action in the gluon sector describe the same thermo¬ 
dynamics as the briefly reviewed NJL model. We con¬ 
clude that in the limiting case of a contact interaction 
both approaches, DS and NJL, give identical results. 
Table 1 contains specific NJL model parameterisations 
taken from Grigorian (2007), reproducing pion and kaon 
masses, and the pion decay constant for different light 
quark constituent masses M^ d . For the purpose of our 
analysis we restrict ourselves to set I and IV with the 
smallest and largest constituent light quark masses in 
vacuum. 


with single-flavor quark and anti-quark spinors (g/,g/) 
and without flavor-mixing terms. In the following para¬ 
graphs, we will not elaborate further on NJL type models 
nor on the techniques to proceed from the Lagrangian to 
describe the thermodynamics of the system. These are 
well described in the literature and for this particular 
model explained in detail, e.g., in Buballa (2005). 

The thermodynamic pressure of a single flavor / is ob¬ 
tained from the grand canonical thermodynamic poten¬ 
tial, 

p = - nf=p FG -^k + ik’ (19) 

with the scalar and vector condensates, <fi and w. Again, 
the free part describes quarks as ideal quasi-particles 
with effective masses, m* = and effective chemical 

potentials, p* = p — Cj. The corresponding NJL model 
gap equations are 

4> = 2G s N c n s (p* f ,m* f ),u> = 2G v N c n v {p},m* f ) . (20) 


4. CHIRAL BAG MODEL 

In this and the following section we introduce a modi¬ 
fied bag model which can be parameterized to reproduce 
NJL model results in very good agreement. We extend 
it afterwards to additionally mimic confinement. 

Fig. 1 shows the dynamical quark masses of u. d , (left) 
and s (right) quarks and the corresponding pressure. 
Chiral symmetry is broken at small and restored at large 
p. The intermediate region is characterized by the ex¬ 
istence of multiple mass gap solutions. The thermody¬ 
namically stable solution maximizes the pressure of the 
system at given p. As indicated by the arrow this re¬ 
sults in a sharp transition from the branch of chirally 
broken (yB) to the branch of chirally restored (yR) so¬ 
lutions. The multiple solutions of our gap equations are 
a feature not restricted to the NJL model. It owes to 
the nonlinearity of QCD’s gap equations in general (cf. 
Chang et al. 2007; Wang et al. 2012; Raya et al. 2013). 
We now fit the pressure of the yR branch by calculating 
only the kinetic pressure of an ideal gas with the bare 
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quark mass to/ and shift it by a constant but flavor de¬ 
pendent value B^. This procedure can be understood as 
a complete transfer of the scalar, regularisation depen¬ 
dent vacuum contributions to the pressure according to 
Eq.(17) into B^. That NJL models in the x'R domain 
can be fitted well by this approach has been shown for 
a flavor blind effective bag constant by Buballa & Oertel 
(1999), and later for flavor dependent bag constants in 
Buballa (2005). 

In analogy to the bag pressure, we name B^ the chi¬ 
ral bag pressure. Evidently, this fit gives good results 
in the %R domain we are interested in. The value of 
B?. is obtained from the condition B£ = P(M°,/x/ = 
0 ) — P(mf, fif = 0), where is the dressed vacuum 
quark mass (cf. Buballa 2005). This prescription has 
first been applied in a DS framework to determine the 
bag constant by Cahill & Roberts (1985). B £ depends on 
vacuum properties only. The critical chemical potential 
at which chiral symmetry is restored in this approx¬ 
imation is defined by the relation Pl in (mf,^Q = B^. 
The model fits yR quark matter according to this pre¬ 
scription for /J > , below jJ x the pressure of the single 

quark flavor and all related quantities can be set to zero 
as this sector is characterized by confined hadrons, con¬ 
sequently the pressure and density of deconfined quarks 
is zero. The existence of a flavor dependent gap Bf is 
a feature not inherent to tdBAG. As we illustrated, it 
originates explicitly from the breaking of chiral symme¬ 
try - an effect which has been deliberately ignored in 
tdBAG as stated by Far hi & Jaffc (1984). Correctly, 
Farhi & Jaffe (1984) refer to the tdBAG bag constant 
as a quantity that mimics confinement and pointed out 
clearly, that the model does not account for chiral sym¬ 
metry breaking. Later in this paper we will point out, 
how these different bag constants are related. The val¬ 
ues of /I x and corresponding chiral bag constants B £ for 
different NJL model parameters are listed in table 2. 


TABLE 1 

NJL PARAMETERISATIONS ACCORDING TO GRIGORIAN (2007). 



M^ d 

[MeV] 

M§ 

[MeV] 

G s A 2 /2 

A 

[MeV] 

m u/d 

[MeV] 

m s 

[MeV] 

I 

330.0 

610.0 

2.17576 

629.540 

5.27697 

171.210 

II 

367.5 

629.6 

2.31825 

602.472 

5.49540 

170.417 

III 

380.0 

636.5 

2.36582 

596.112 

5.54297 

169.559 

IV 

400.0 

647.7 

2.44178 

587.922 

5.58218 

167.771 


5. VECTOR INTERACTION IN THE BAG MODEL 

Before we address this issue in detail we extend the 
chiral bag model to account for the repulsive vector in¬ 
teraction. This aspect is crucial if one aims to describe 
stable and massive hybrid neutron star configurations 
with quark matter cores, cf. Klahn et al. (2013), in agree¬ 
ment with the recent observations of pulsars with two 
solar masses by Demorest ct al. (2010) (PSR J1614-2230 
with 1.97 ±0.04 M 0 ) and Antoniadis et al. (2013) (PSR 
J0348-0432 with 2.01 ± 0.04 M 0 ). 



g [MeV] 

Fig. 2. — (color online) Impact of vector interactions on the 
stiffness of the EoS and the corresponding shift of the critical po¬ 
tential ji x which defines the chiral phase transition. Thin lines: 
vBag. Bold lines: NJL model. Parameters: set IV, table 1. 


The assumption that the dressed quark mass in the 
X'R region is effectively equal to the bare quark mass to 
and the pressure of the system evolves mainly due to 
the corresponding kinetic contribution while all effects 
of the scalar condensate are hidden in the constant off¬ 
set Bf, frees us from solving the mass gap Eq. (9), the 
only equation which had to be solved self consistently. 
The actual chemical potential /i is easily obtained from 
Eq. (10) which ensures that the vector gap Eq. (5) is 
solved correctly. It requires to perform only a single in¬ 
tegration in order to compute the one-particle number 
density n/(/rJ) = n„(/x)., to/). All further calculations 
are simple rescaling operations involving free Fermi-gas 


expressions in terms of the effective chemical potential 
/i*. The rescaling results from the vector-gap equation 
Eq.(20) and affects the thermodynamics in terms of pres- 

sure and energy density, 


Af = + K v n f (fj,* f ) , 

(22) 

Pffaf) = PfgM) + §^ 2 /(/4) - K , 

(23) 


(24) 

3 

II 

3 

>-* 

~C 

(25) 


These equations (except for a few remaining consid¬ 
erations regarding deconfinement in sec.6) define the 
novel vBAG model and are sufficient to reproduce NJL 
model results with finite values for the vector coupling. 
For convenience we introduced a new coupling constant 
K v = 2G V . In Fig. 2 we reproduce NJL model results 
for different values of the vector coupling constant Gy 
at a fixed value for the scalar coupling constant Gs and 
compare them to the bag model approximation. As the 
scalar interaction is unchanged the vacuum offset Bf is 
the same for all shown scenarios. The critical chemi¬ 
cal potential shifts to higher values because of the 
increasing vector coupling which increases the chemical 
potential according to Eq. (22). This effect dominates 
in comparison to the additional contribution to the pres¬ 
sure (2 nd term in Eq. (23)) which just by itself would 
increase P at given /i* and therefore shift towards 
smaller values. 
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Fig. 3. — (color online) Symmetric matter phase transition from 
nuclear to (u,d)-flavor quark matter using the nuclear model EoS 
TMl(blue line), from Shen et al. (1998), NJL model for 2-flavor 
quark matter using the same parameterization as in Fig. 1 (black 
circles) and the bag model parameterization of the NJL result (red 
lines). Applying the standard Maxwell construction the crossings 
with TM1 are thought to mimic the deconfinement phase transi¬ 
tion. Red dashed: By default, the so defined deconfinement chemi¬ 
cal potential [idc is larger than ji x . Red dash-dotted: The pressure 
is shifted by a positive value = Ptmi (l l x ) to match \i,\ r = fix 
when the pressure of quark and nuclear matter equal. This mod¬ 
els simultaneous deconfinement and chiral symmetry restoration. 
Parameters: set IV, table 1. 

6. CONFINEMENT TRANSITION 

The model as introduced so far is suited for deconfined 
quarks but as typical for NJL models it cannot describe 
truly confined matter 2 . Already accounting for bound 
states is a demanding problem that requires analyses of 
two- and three-particle in-medium correlations. Only few 
studies exist which aim to describe dense hadronic matter 
explicitly in terms of quark matter degrees of freedom, 
e.g., by Wang et al. (2011) and Blaschkc ct al. (2014). 
For phenomenological purposes deconfinement is there¬ 
fore usually modeled via a first order Maxwell transition 
from a hadronic to a quark matter EoS which both are 
not modeled within the same framework. The critical 
chemical potential Hdc which defines the deconfinement 
transition is obtained by applying the Maxwell-condition 
Ph(h< 1 c) = Pq{h<Ic)- We illustrate this in Fig. 3 for 
isospin-symmetric matter at zero temperature , where the 
nuclear EoS is the relativistic mean-field model TM1 of 
Sugahara & Toki (1994) and Shen ct al. (1998). A TM1 
based EoS by Hcmpel & Schaffner-Bielich (2010) can be 
found online in the CompOSE-database 3 . CompOSE 
is explained in detail in Typel et al. (2013). The follow¬ 
ing discussion does not rely on TM1 as a specific choice; 
our arguments hold in general and the same qualitative 
results would be obtained for any different nuclear EoS. 

For the quark matter branch we neglect the strange 
flavor and consider a purely chiral bag model for iso-spin 
symmetric two-flavor quark matter ( n u = rid)- With 
the chosen parameterization we obtain the two-flavor 
chiral bag constant B™ ,d = B™ + B d = (181.6 MeV) 4 , 
which results in a chiral critical chemical potential of 
p x = 1148 MeV. The Maxwell-condition without further 
modifications to the model is then fulfilled at a decon- 

2 However, confinement understood as the absence of quark pro¬ 
duction thresholds can be mimicked even in the non renormalizable 
NJL model (cf. Gutierrcz-Guerrero et al. 2010). 

3 http://compose.obspm.fr/spip.php?article29 


fmement critical chemical potential of Hdc = 1651 MeV. 
The intermediate region spans an interval of about 500 
MeV in which the physical interpretation is ambiguous. 
Our main concern is that the QM model predicts the 
restoration of chiral symmetry which should be reflected 
by the nuclear model EoS which is of course not the 
case, as TM1 is an independent model. There are sev¬ 
eral possibilities to deal with this situation. The first 
would be to accept it as a poor description of hypo¬ 
thetic quarkyonic matter where the chiral symmetry of 
light quarks is restored in a gas of still confined nucleons 
(cf. McLerran & Pisarski 2007; Glozman & Wagenbrunn 
2008, and references therein). Microscopic studies indi¬ 
cate, that the transitions related to chiral restoration and 
deconfinement are closely related and can occur at the 
same or very similar chemical potential (cf. 3ender et al. 
1996, 1998; Blaschke et al. 1998; Qin et al. 2011). In or¬ 
der to illustrate the setup of our fully defined model we 
prefer such a situation and assume that chiral restora¬ 
tion and deconfinement occur at the same critical chem¬ 
ical potential, hence p x = pdc- We model a hybrid 
EoS which fulfills this criterion and introduce the addi¬ 
tional deconfinement bag constant, Bd c , which is added 
to the overall pressure of our quark matter model such 
that Ph(P' X / dc) =: PQ^x/dc)- The same strategy has 
been suggested in Pagliara & Schaffner-Bielich (2008) as 
a sound alternative to just perform a standard Maxwell 
phase transition from a nuclear model EoS to a vacuum 
renormalized NJL quark matter EoS. While we consider 
this as a likely scenario we point out, that a rather 
quarkyonic behavior with p x ^ pdc is not ruled out. 
A strength of our model is the independent treatment 
of the deconfinement transition and chiral restoration, 
where the latter is defined only by the chiral bag con¬ 
stants and the relation -P/ irl (m/, P x ) = B^, whereas the 
deconfinement critical potential depends on all chiral bag 
constants B{ and the nuclear model EoS. The total pres¬ 
sure in our model is P = ^2f(P^ m +K v nj/2)—B e ff with 

B eff=J2 B X~ Bd - ( 2 6 ) 

/ 

We emphasize, that a single quark flavor can contribute 
to the sum of chiral bag constants only if chiral symmetry 
for this flavor is restored; as explained in sec. 4. This is of 
particular importance for the s-quark if /i s > ^u,d Notg, 
that this prescription - subtracting an effective bag con¬ 
stant from the sum of all partial kinetic pressures - corre¬ 
sponds to the standard tdBAG prescription. In addition 
to tdBAG, our model identifies different contributions to 
the effective bag constant. This is not a minor technical 
detail. One has to keep in mind, that the original tdBAG 
explicitly omits the scalar (and vector) interaction and 
attributes the bag constant fully to a confining back¬ 
ground field with positive energy and hence negative, 
confining pressure. In our prescription, the positive value 
of the bag constant results only from the restoration of 
chiral symmetry, while confinement/deconfinement, al¬ 
though introduced merely phenomenologically, reduces 
this value. From our perspective, this makes perfect 
sense if one recalls, that the most naive perception of 
confinement is the binding of quarks in the chiral broken 
phase. Thus, our deconfinement bag constant is easily in- 
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terpreted as the actual binding energy of confined quarks 
which effectively reduces the free energy of the system. 

7. NEUTRON STARS WITH QM CORE 

In this section we discuss how vBAG clears the long 
lasting perception originating from tdBAG, that quark 
matter tends to be too soft towards higher densities in 
order to support the idea of QM in compact stars. 



Fig. 4. — (color online) vBAG EoS pressure vs. energy den¬ 
sity for neutron star matter (left panel) and corresponding mass- 

radius relations (right panel) for ^ = 137.6 MeV and 

= 224.1 MeV (set I, table 2), comparing different values of 
Gy. Values of the deconfinement bag constants for the different 
Gy are listed in the left panel. NJL parameters: set I, table 1. 

In the left panel of Fig. 4 we illustrate the phase tran¬ 
sition in /3-equilibrated neutron star matter for selected 
chiral bag constants, B 1 ^^ ^ = 137.6 MeV and B\[x = 
224.1 MeV, and varying vector interaction strengths 
( G v = (0,1/4,1/2)G S ). Bdc has been adjusted such that 
= l^dc- /J-x varies for different values of the vector cou¬ 
pling G„, as illustrated by the inset of Fig. 2, and there¬ 
fore the value of Bd c depends on G v as well. Variations 
of the s-quark onset density are almost negligible. 

Due to the large vacuum mass of the s-quark, the phase 
transition from TM1 to two-flavor (it, d ) quark matter 
takes place at lower density, followed by the transition 
to three-flavor (u, d, s ) matter at higher density. This is 
the behavior one expects from NJL-type models without 
flavor coupling channels. It is not accounted for by td¬ 
BAG which ignores DySB and consequently predicts a 
transition from nuclear to three-flavored strange matter. 
In contrast, vBag describes a sequential transition from 
nuclear to two-flavor, then to three-flavor quark matter, 
where the three flavor branch of the EoS is found only for 
unstable NS configurations, see the right panel of Fig. 4. 

Note that for the canonical vector coupling of Gy = 
1/2 Gs vBAG nearly reproduces the nuclear model 
EoS TM1. This “masquerade” has been discussed by 
Alford et al. (2005). It results in a very similar mass- 
radius relation for pure neutron and hybrid stars for this 
parameter set (right panel of Fig. 4) . The different 
onset-densities for two-flavor matter at different values 
of G v are reflected in different neutron star mass-radius 
relations. For G v = 0 the transition in general results in 
unstable configurations for the parameters explored here. 
Hence, the vector interaction is essential to stabilize the 
compact stellar object. Larger values of B^ d associated 
with larger quark masses result in higher critical densities 
for the phase transition but qualitatively reproduce the 


above discussed features as long as the transition density 
does not reach values where already the purely nuclear 
NS configurations render unstable. 

8. STABILITY ANALYSIS OF STRANGE MATTER 

The long-standing hypothesis of absolutely stable 
strange quark matter as the ground state of strongly in¬ 
teracting matter, introduced by Witten (1984), is sup¬ 
ported by the standard tdBAG of Farhi & Jaffe (1984). 
As stated before, tdBAG does not account for DySB nor 
the vector interaction channel. NJL model studies imply 
that taking DySB into account absolutely stable strange 
matter can be ruled out, see e.g. Buballa & Oertel 
(1999). Our model is a phenomenological hybrid of NJL 
and tdBAG and, in terms of the defining parameters, 
namely BBd c and K v , could in principle be used to ar¬ 
gue against or in favor of absolutely stable strange matter 
if one decides to choose the corresponding values freely. 

Briefly described, the concept of absolutely stable 
strange matter is based on the idea, that two flavor 
matter at zero pressure has a larger energy per par¬ 
ticle than the most stable nuclei (we chose a value of 
931 MeV) because otherwise these would decay into the 
constituent u- and d- quarks. Adding s-quarks to the 
mixture, however, lowers the energy per particle. Hence 
strange matter could be more stable than nuclear matter 
(cf. Farhi & Jaffe 1984). This leaves a small window of 
possible (small) bag constants which support this idea. 
The main reason why chiral models do not confirm this 
idea is, that the chiral bag constant itself is flavor depen¬ 
dent due to DySB and does not differ from the 2-flavor 
bag constant due to the large value of the dressed s-quark 
mass and the resulting high densities where free s-quarks 
can appear. More precisely, tdBAG predicts absolutely 
stable strange matter in a domain where chiral models 
predict the chiral symmetry for the s-quark to be broken 
and to be restored for the light quark. This is the defini¬ 
tion of two flavor quark matter where s-quarks do not ex¬ 
ist. This situation holds in the vertical gray shaded band 
(DySB) in Fig. 5. In the domain where the ^-symmetry 
of the s-quark is restored the energy per particle of the 
2-flavor phase is already high and would be increased fur¬ 
ther by the chiral phase transition. This is illustrated by 
the red solid line in Fig. 5 with B ® where /n® is 

indicated by the black vertical dashed line near 2 GeV. 

In general, the chiral bag constant as obtained from 
NJL model analyses has a large value of > 160 MeV 
already in the two flavor case, a value which exceeds the 
one predicted for stable strange matter. Therefore, we 
use the additional free parameter Bd c to reduce the ef¬ 
fective bag constant and the energy per particle. Here 
we varied Bd c = (0,1,2,3) times B^ d ■ In order to repro¬ 
duce the assumptions of tdBAG we set = 0 (black 
lines in Fig. 5). It lowers the energy, however, only sig¬ 
nificantly at [i < hence a region governed by DySB 
of the light flavors where matter is still confined and a 
model for deconfined matter does not apply. Note that 
the standard tdBAG would still predict a scenario with 
stable strange matter as it does not account for DySB. 
vBAG predicts effectively the same bag constant but ex¬ 
cludes the existence of free quarks from this region, since 
the chemical potentials are below the critical chemical 
potential for DySB, ^ < /r x . 
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Fig. 5. — (color online) /3-equilibrated strange quark matter: 
energy for different values of Bdc at T = 0, compared with the 
nuclear matter minimum (horizontal shaded band). The chiral 
bag constant B x has been chosen from the parameter set I, table 1 
with a dressed light quark mass (in vacuum) of 330 MeV to provide 
a number close to the lower limit of B x . red: finite B hence 
there are two sequential chiral transitions (nucl.-f^-2f-f^3f). black: 
B J = 0 results in only one chiral transition (nucl.-f^3f). Note that 
the vertical shaded band marks the domain of D^SB which tdBAG 
cannot describe. See text for more details. 

The last available parameter is K v . From Eq. (24) 
follows, that the vector interaction adds to the energy 
density and therefore makes strange matter less stable 
at given density. Moreover, as discussed earlier, it shifts 
the chiral transition to higher chemical potentials which 
further disfavors stable strange matter. 

Back to our analysis, we conclude that the idea of sta¬ 
ble strange matter is not supported for a NJL based pa¬ 
rameterization of the chiral bag constants, even after in¬ 
troducing a deconfinement bag constant, which leads to 
effective bag constants similar to the original tdBAG. 

How general is this result? The DS study of 
Nickel et al. (2006) implies that the color flavor locked 
phase (where all three quark flavors exist) is favored over 
the 2SC phase which consists of the two light quarks only. 
Would this make absolutely stable strange matter more 
likely? We claim that not. The crucial question regard¬ 
ing the stability of strange matter is the onset of chiral 
symmetry restoration in terms of the chemical potential 
for all three flavors. This is, because absolutely stable 
strange matter is predicted in tdBAG as a chiral sym¬ 
metric phase. As soon as the light quarks attain mass 
due to DySB the corresponding baryon energy per par¬ 
ticle has to be either higher (no confinement) or equal 
(confinement) to that of nuclear matter. If it would be 
lower, the argumentation of Far hi & Jaffe (1984) holds: 
nucleons would simply decay into their constituent light 
quarks. Our simple NJL model predicts DySB for the 
light quarks at critical beyond values which support 
the absolutely stable strange matter hypothesis. In con¬ 
trast, Nickel et al. (2006) showed that the current mass 
of the s-quark and therefore p,® reduce due to color-flavor 
locking. However, the color-flavor locking at the same 
time increases /z“’ d , in fact they coincide, /z® = /z“ ,d . We 
argue, that this increase of the light quark critical chem¬ 
ical potential for DySB acts even stronger against the 
formation of absolutely stable strange matter. 

We point out again, that the crucial mechanism which 
prevents the existence of absolutely stable strange matter 
is DySB in the light quark sector, an explicitely feature 


inherent to non-perturbative QCD. Therefore, any anal¬ 
yses of this hypothesis that bases on a perturbative ap¬ 
proach to QCD, e.g., Kurkela et al. (2010), is not suited 
to derive reliable qualitative conclusions. 

9. SUMMARY 

Quark matter in compact stellar objects, such as 
(proto)neutron stars, has been studied widely based 
on the simple and powerful tdBAG model. With re¬ 
cent advances in mass-determinations of massive com¬ 
pact objects it has become increasingly difficult to fulfill 
the ~ 2 Mg maximum mass constraints with tdBAG 
quark matter EoS. It becomes particularly difficult if the 
hadron-quark phase transition yields large latent heat, 
i.e. large jump between hadron and quark EoS. In this 
article we define the novel vBAG model by introducing 
vector interaction. Since vector interactions stiffen the 
EoS with increasing density this provides an undisputed 
cure to the aforementioned problem without raising the 
technical difficulties in computing the EoS. 

We derived vBAG from the currently most sophisti¬ 
cated approaches to dense quark matter. We demon¬ 
strate how a simple contact interaction on the gluon 
propagator level allows us to reproduce NJL and Bag 
models within the DS formalism in medium (i.e. at fi¬ 
nite chemical potential). Moreover, we illustrated how 
the vector interaction channel is a crucial input for the 
investigation of hybrid, quark matter - neutron stars. 
Confinement is mimicked by a confinement bag constant 
Bdc■ It reduces the sum of all flavored chiral bag con¬ 
stants. Together, they reproduce effective bag constants 
typical for tdBAG (cf. Farhi & Jaffe 1984). We revisited 
the absolutely stable strange matter hypothesis and con¬ 
clude that it cannot be easily discussed without a thor¬ 
ough understanding of DySB in particular of the light 
quarks. Based on our extended bag model, which im¬ 
plements DySB in a coarse fashion we confirm previous 
NJL model studies that rule out the existence of stable 
strange matter and that tdBAG’s prediction of stable 
strange matter results directly from the suppression of 
D\SB. Confinement in terms of Bd c is not likely to be 
an effect which would change this result. 

The model EoS vBAG has been deliberately developed 
to be simple and easy applicable while catching the main 
(minimal) features we expect from QCD in dense matter, 
namely DySB and to some extent deconfinement. We 
consider it an important improvement of the standard 
bag model as it tries to clearly separate these effects. Al¬ 
though the number of parameters is higher, each of them 
has a clear physical interpretation. Note, that QCD in 
medium is not plainly solvable, which was reminiscent 
when we motivated the model from the DS perspective. 
We think that vBAG is a practical tool for modelers who 
wish to account for QCD degrees of freedom in complex 
dense systems - in particular for applications in astro¬ 
physics studies. 
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TABLE 2 

Single flavor, effective two-flavor chiral bag constants and p x / dc for the parameterisations of table 1 . 



chiral bag model parameters 


Phase Transition TM1 —> 2f QM (symmetric) 


i 

pu 4 

r BAG 

(E-Pbag) 4 

1 

ps 4 

r BAG 

( E ^sag) 4 

u/d 

the 

/4 

Vx 

Mdc (Maxwell) 

pTMlfc) 

(x-dc) 


[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

[MeV] 

I 

137.6 

163.6 

224.1 

238.5 

343 

594 

1029 

1458 

120.7 

149.8 

II 

145.8 

173.4 

221.9 

240.2 

365 

591 

1094 

1569 

141.4 

149.8 

III 

148.5 

176.6 

221.7 

241.3 

371 

590 

1114 

1600 

147.1 

149.9 

IV 

152.7 

181.6 

221.7 

243.3 

383 

590 

1148 

1651 

155.3 

150.0 
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